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Abstract 

The d - dimensional n - colour lattice h d with modular sublattices are studied, when 
the only one crystallographic type of sublattices does exist and the only one of the colours 
occupies a sublattice, which is still invariant under fc-fold rotation Such kind of colouring 
always preserves an equal fractions of the colours composed h d . The n - colour lattice with 
modular sublattices allow to exist the crystallographic rotations Ck for k = p r , r > 1 and 
n < p, where p is a prime number. 



I. INTRODUCTION 



One of the fascinating results in the d - dimensional (dD) crystallography is the Hermann 
theorem on the crystallographic restriction (CR): in dD - lattice L d there are allowed the 
crystallographic rotations Ck of all k which satisfy d>^f(k), where ty(k) is the Euler totient 
function. Due to Schlaffi [Q, a lattice L M+1 of odd dimension always possesses a rotation 
axis Ck if a lattice L 2d possesses this axis, that reflects an even values of the totient function 
#(Jfe). 

Recent progress in the problem of effective conductivity in 2D three - component regular 
composites || poses a question about CR in a plane tiling perfectly coloured with n different 
colours of equal fractions. This question is concerned with another Hermann theorem 0] 
about material tensors of r-th rank operating in media with an inner symmetry that includes 
a rotation axis of the highest order I. namely, the existence of the rotation axis C\ of order 
I > r makes the r-th rank tensor isotropic. There are two other sorts of colour lattices 
which are of high importance in physics of condensed matter. These are the lattices with 
colour edges, which appear in the bond percolation problem [|[ and the lattices with colour 
nodes which are related to the site percolation problem in the theory of phase transition ||. 

The last kind of colour lattices is most simplest from crystallographic standpoint. Re- 
cently the group - theoretical methods were successfully used here to find the different 
possibilities to colour the lattice sites 0, ||. We have in mind the coincidence problem 
for the colouring a lattice possessed a k - fold rotation with finitely many pairwise distinct 
colours in such a way that the colour version is still crystallographic: one of the colours 
occupies a sublattice which is still invariant under k - fold rotation, while the other colours 
label the cosets. In other words a problem was reduced to the question: how many sublat- 
tices with k - fold rotation of h d of a given index n do exist. In 2D lattices the classification 
of such coincidence is the same as the classification of the colour lattices with rotational 
symmetry ||. For such kind of colouring the obtained results 0, || can be interpret as a 
generalization of Hermann theorem for colour lattices, e.g. in the 2D n - collour lattices the 
following crystallographic rotations Ck are allowed 

C 4 _> n = 1, 2, 4, 5, 8, 9, 10, 13, 16, ... , C 6 -> n = 1, 3, 4, 7, 9, 12, 13, 16, ... . 



For indecomposable point groups in higher dimensions no nontrivial orthogonal transfor- 
mation commutes with all symmetries so there is no longer this equivalence. One can see 
however that even in the 2D - case the colouring rules defined in 0, || admit both the 
lattices with a permutation invariance of different sublattices as well as the lattices where 
this property is absent (see Fig. [I]). Moreover such colouring admits a non-equal fractions 
of the colours composed the entire lattice. 

A new algebraic approach [|K]] to the effective isotropic conductivity problem of 2D 
infinite non - random composite has shown the existence of self - dual algebraic functions 
A n (#i) which are invariant upon the action of full permutation group S n . The regular 
arrangement of n colour nodes at the plane, which are corresponding to distinct constituents, 
is restricted by growing n: the existence of rotation symmetry with Ck, k > 3 preserving the 
permutation invariance of the equal fraction constituents disappear when n > 4. By other 
hand the mentioned properties can be restored if we admit the quasi-crystalline arrangement 
of the nodes at the plane via their projection from the high-dimensional colour lattice, which 
possessed a high-fold rotation axis. 

This motivates us to study such kind of colouring where the only one crystallographic 
type of sublattices does exist: any two sublattices L^(n) and lA(n), which belong to the 
distinct colours q± and q<i respectively, are isomorphic IJzAn) «-> U^ 2 (n) via translational 
motion by a distance and in a direction equal to one of the basic vectors. An isomorphism 
preserves the equal fractions of the colours in the unit cell. Such sublattices will be called 
modular due to their relation to the modular equations. 

The colour lattices, which we are going to study in the present paper, give rise to another 
version of Hermann theorem: what rotational axes are survived after colouring the nodes 
of the lattice h d . In the case d = 2 an answer was found in [IT| where an existence of the 
crystallographic rotations with k = 2,3,4 for 2 - colour lattices and k = 3 for 3 - colour 
lattices was shown. This means that the isotropic tensor of second rank, like conductivity 
a, exists in the plane media regularly composed of not more than 3 sort of nodes. 

The approach which we exploit here is quite simple: starting from the admitted axes Ck 
we will check which of them are survived after colouring the lattice in n colours. 



II. TWO DIMENSIONAL LATTICES 



Before going to CR in the colour lattices we recall its derivation for 2D unicolour (n — 1) 
lattice. A plane lattice L 2 is arranged as a set of periodic sequences of nodes defined by 
vectors d 2 G L 2 

d 2 = m\&i + m 2 e 2 , (e 1 , e 2 ) = cosip, ^^0,ir, 

where m^mj £ and e!,e 2 are two noncollinear unit vectors, called lattice basis. Such 
set contains infinite number of points. Any translational motion of this lattice by a distance 
and in a direction equal to one of the basic vectors superimposes the points of the set upon 
their original positions, so that the new state of the set is indistinguishable from the old 
one. 

The rotational part R of every crystallographic symmetry operation has a representation 
through a unimodular 2x2 matrix with integral entries. It is simply to establish this matrix 
for counter-clockwise isometric rotation R of the order k = 2tc / (fik , (j^kj — h '■ 

( 



Rk 



sin 



\ 



sin(^ fc - <p k ) -sin0 fc 
sin0 fc sin(V> fc + fc ) 



, det Rk 



(1) 



where J 2 is a 2 x 2 identity matrix. The admitted values of ipk and <f>k are exhausted by the 
following list 



Ci : 0i = 27T, R 1 = I , C 2 : 02 = tt , R 2 — —I , 

-1 , k = 3 



2tt 



Cfc : <f>k — i>k — -Rfc 



^0 -1^ 



y 1 Ofc y 



, fc = 4 (2) 

1 , k = 6 

We will define now a colour plane lattice L 2 (n) as a union of n modular sublattices L 2 (n), g = 
0, 1, n — 1, each of them has a separate colour. The construction of q-ih sublattice L 2 (n) 
is the following 



h 2 q (n) = |J d 2 ^, 0<g<ra-l, 

d 2 9 ' n = m^ei + m^e 2 , mS ?) + = g mod (n) 



(3) 



and all cyclic rotations [R k ) preserve the invariance of sublattice L 2 (n) 



R k j d 2 q ' n e h 2 q (n) , 0<t< fc-1. (4) 

Here M q>n denotes a set of integer solutions , m^' of the system of A; linear modular 
equations (|3], [|). Such definition of L 2 (n) has one important corollary: if L 2 (n) is invariant 
upon the isometric rotation R k then there exist [n — 1) - colour lattice L 2 (n — 1) which is 
also invariant under rotation R k . Indeed, let us delete one modular sublattice Lf(n) from 
the entire L 2 (n). This does not affect the existence of the other n — 1 relations (|3|, [|) and 
therefore the h 2 (n — 1) - lattice remains to be invariant. In fact, this is true for any number 
of deleted sublattices. Thus, we conclude that if the L 2 - lattice could be coloured in n 
colours, then it also could be coloured in any number smaller than n. We denote such 
maximal number of colours in the plane lattice with N%. 

Let us consider now a plane lattice L 2 spanned on the basic vectors ei,e 2 



ei 



V / 



e 2 



V) 



and possessed a rotation axis Ck- We mark the lattice in n colours in accordance with 



(|j). Acting with an isometric rotation Rk (0) onto the entire colour lattice L 2 



n 



Rk ei = e 2 , R k e 2 = - e x + a k e 2 



we act by the same on all sublattices L 2 (n) 



Rk d 2 



q,n 



(q) 



—777.2 e l + ( m l + a fc m 2 ) e 2 



(5) 



Let us consider the corresponding system of two modular equations t = 0, 1 



nil + m 2 = Q mod (n) , q — 0, n — 1 
+ (a^ — l)rn 2 = ? mod (n) . 



(6) 
(7) 



Both equations (H), (^) are coincident when the equality (2 — a k ) = mod (n) is valid. This 
leads to the conclusion 



N 2 k = 2 - a k 



N% = 3 , iV 2 = 2 , JVf = 1 



(8) 



in accordance with [|ll[| . It is simply to convince that the other k — 2 modular equations 
do not contradict with ([8]). 

At Figure we present two different colourings of the plane lattice with modular sub- 
lattices which preserve the C3 rotation axis. The similar treatment for two - fold axis C 2 is 
trivial 

R 2 ei = - ei , R 2 d 2 9 >" = -d 2 *>" , 

that yields 

mf ] + m { 2 q) = ±q mod (n) — ► 2 (mf ] + m 2 q) ) = mod (n) — ► iV 2 2 = 2 . (9) 

In 3D lattice the numbers of colours N% which preserve the modular invariance (0) remain 
the same like in 2D lattice. 



III. HIGH DIMENSIONAL LATTICES AND UNIMODULAR MATRIX 

REPRESENTATIONS 



In this Section we consider CR in high - dimensional colour lattices h d (n) with modular 
sublattices. We define such lattice in a way similar to @: it is a union of n modular 
sublattices h^(n),q = 0, 1, ...,n — 1, each of them has a separate colour. The construction 



of g-th sublattice IL^(n) is the following 



u 



d g J n , 0<q<n-l 



m\ q) £M q , n 



dr = 5>* (9)e i> E 

i=l i=l 

where the basic vectors read 

(A 



m 



(</) 



q mod (n) , < q < n — 1 , 



ei 



, e 2 



v / 



, e d 



\ / 







(10) 



A straightforward way to built out the linear modular equations using the matrix representa- 
tions for every axial rotation Cf. in h d is cumbersome due to rapid growth of both a number 
k of equations and a dimension of irreducible representations. To avoid such algebra 
we will make use of a special property of axial rotation C p , where p is a prime number. This 
is a unique rotation, all the powers of which have irreducible representations. Such sort of 
rotations arranges a natural basis to utilize all other rotations with non - prime k. 



A. Rotation symmetry C}~, k = p 



This symmetry operation has an irreducible representation R„ through a unimodular 



companion matrix [12J det R p = 1, dimR p = p — 1 



R,, 



... 

1 ... -1 
1 ... -1 



... 1 -1 
... 1 -1 



FL 



1 , if j = i- 1 , 2 < i < p- 1 
-1 , if j =p- 1 , 



, otherwise , 



where f R 



= I n , diml n = p — 1 and 1^ is an identity matrix. In other words R p is a 
generator of Abelian group C p over the integer numbers Z. All the powers of R p are also 
irreducible matrices 



R, 



1, ]fj=p + i — t, 1 < i < t - 1 
-1 , if j = p - 1 , 

1 , Hj = i-t, t+l <i <p - 1 
, otherwise . 



(12) 



E.g. 



/ 



Re 



-1 \ 
10 0-1 
10-1 
1 -1 J 



( 



Rr 



0-1 1 \ 
0-10 
10-10 
1 -1 Oy 



/ 



R, 



110 
•10 10 
•10 1 
•10 



\ 



(13) 



Let us present its action on the basic vectors ei 



or more generally 



r>— 1 

= -J2 e i > 


-Rp e j = e j+i , j = i,-,p-2 • 


(14) 


(■^p) e p-t-j 


= e P-j > J = l,...,p-t - 1 , 




^i? p j e p _ t = 


p-i 

-J^e; , t = l,...,p-l , 

i=l 




1 R p ) e p _ t +j 


= ej , j = 1, — 1 . 


(15) 



Then the action on a generic vector dp_i is 

j=t+l j=l j=l 

i.e. 

J=2 j j'=3 i=l 

The modular equations are the following 

p-i 

= g mod (n) , g = 0, n — 1 , (16) 

j=i 

p-i 

— p rnf"' = q mod (n) , i = 1, — 1 . 
They are coincident if p = mod (n). Hence it follows N^ 1 = p. 



B. Rotation symmetry Ck, k = P1P2, Pi / P2 • 

The case k = P1P2, Pi 7^ P2 presents another situation, when the modular equations do 
not have any solutions. First we show this in the trivial case k = 2p, p 7^ 2. The irreducible 

representation R 2p = —R p for axial rotation C 2p leads finally to the modular equations 

p-i 

rnf = ±q mod (n) , g = 0,...,n-l, (17) 

j'=i 
p-i 

^ra|''-pm-'' = ±5 mod (n) , i = l,...,p— 1. 

j'=i 



or, more briefly, 



p-i 



('/) 



p = mod (n) 



which do not have solutions for any modulus n: N^ p l = 1 . 

The generic case p\ 7^ P2 leads to two equivalent irreducible representations through a 
{pi — \ ){jp2 — I) — dim unimodular matrices A piP2 = R P1 <g> i? P2 and -B P1P2 = i? P2 <8> R P1 , where 
(g> denotes a tensorial product. It is convenient to present they as block - matrices with 
companion matrix entries 



/ 



A 



pipy 



0"7Ti 0"7Ti 0"7Ti ••• 0"7Ti 0"7Ti 



\ 



Pi 



■^pi ••• O^i Otti 



-it 



Pi 



0?ri i^pi Ovri ••• 7ri 0^ i^pi 



V 



O^i ••• i^pi Otti i^pi 



Ojri O^i O^i ••• O^i i^pi i^pi 



0^2 ^71-2 0^2 ••• 0-7T2 i^: 



) Rplp2 



l-p 2 



Rp2 0^2 0^2 ••• 0vT9 



VT2 u vr2 



-it 



P2 



0-7T2 i^P2 ^7T2 * * * ^7T2 ^7T2 i^P2 



0-7T2 ^7T2 ^7T2 ' ' * i^P2 ^7T2 i^P2 



0-7T2 0^2 0^2 "• 0^2 i^P2 i^P2 J 



where R Pl , R P2 can be extracted from (JTTJ) . The zero - matrices 7ri , n2 have dim 7ri = 
Pi — 1, dim 0^ =p 2 — 1. The successive powers of A piP2 and -B piP2 are the diagonal block - 
matrices 

P2-1 



.4 



P1P2 



B 



P1P2 



P2^2 



Pl*l 



diag 



diag 



R 



P2*2 



1 R Pl 



P2*2 



1 i^pi 



P2*2 



Pl-1 



7? 



P2 



Pl*l 



' I -^P2 



Pl*l 



I i^P2 



Pl*l 



*2 = l,-,Pl 



ti = l,...,pa • 



(18) 



We have one more simplification in the last formula? which comes due to gcd(pi,p2) = 1 and 



it 



pi 



pi 



ijri > ( i^ 



l P2 



P2 



I n2 . Indeed, when a variable t 2 runs through the positive integers 

/ ^ \ P2*2 / ^ \ *2 

1, pi then the matrix representation ( i? Pl j runs actually through the matrices ( R P1 
and two sets of matrices coincide, e.g. for p x = 5, p 2 = 7 we have 



R 5 ) 7 = (r 5 )\ (r 5 ) U = (r 5 )\ (r 5 ) 21 = (rX , (r 5 ) 28 = (r 5 ) 3 , (r 5 ) 35 = (rX 



The same is valid for the running variable t% and matrices ( R 



'P2 



Pi*i 



i.e. 



R 



pi 



P2*2 



it 



PI 



i? 



P2 



Pltl 



R 



P2 



(19) 



Thus the diagonal matrices fll8|) can be rewritten as following 



A 



\v2t2 



PVP2 



j = diag 



*2 



Rpi ) ) ( Rpi ) )■■•;( Rpi 



pi-i 



t 2 



B 



P\P2 



Pltl 



diag 



Rp2 ) ) ( Rp2 ) )•••)( Rp2 



\ P2*2 



t 2 = 1,...,P! 



, tx = 1, ...,p 2 



(20) 



Let us consider the action \A PlP2 j on the basic vectors ej. It can be decomposed in p 2 — 1 
parts s = 0, . . . ,p2 — 2 



*2 

-Rpi ) e s(pi-l)+pi-t 2 -j = e s(pi-l)+pi-j , h = 1, •••,£>! - 1 



R„ 



Pi j e s(pi-l)+ 



(s+l)(pi-l) 



Pl-*2 



i=l+s(pi-l) 



Rrs 



pi j e s(pi-l)+pi-t 2 +j — fs(pi-l)+j 



j = 1, -,Pi - h - 1 



(21) 



The corresponding modular equations are 

(P2-l)(pi-l) 

= g mod (n) , g = 0, — 1 



(P2-1)(P1-1) 



P2-2 



mf -Pi ^2 m { ^ (pi _ 1} = q mod (n) , i = 1, - 1 . 



A similar consideration concerned with the action ( B PlP2 ) leads to the modular equations 



P2<2 



(22) 



(p 2 -l)(pi-l) 

E 

i=i 

(P2-1)(P1-1) 

E 

i=i 



(9) 



q mod {n) , g = 0, n — 1 

Pi-2 



(9) V" 

s=0 



ra- 



il) 



g mod (n) , i = 1, ...,P2 — 1 • 



Both equations fl22|) and (|3|) give a system 



(23) 



P2~ 2 Pi" 2 I 7 = 1 m - 1 

s=0 s=0 I j = 1, ...,p 2 - 1 . 



(24) 



which do not have solutions for any modulus n: N P ^ P2 1 ^ P2 ^ = 1 . In Appendix [A] we 



present two irreducible representations Am, B^ for C15 rotation and give a straightforward 
solution for corresponding modular equations. 



C. Rotation symmetry Ck, k = p r , r > 2. 



This symmetry operation has an irreducible representation R p r through a unimodular 
companion matrix |T2j det R p r = 1, dim R p r = p r ' 1 (p — 1) 

' ... -l\ 
1 ... 



Rr> 



1 ... 

... -1 

... -1 

... 1 

... 1 



'•J 



1, if j = i- l, 2<i< p r ~ l {p - 1), 
■1 , if j = p r ~ 1 (p — 1) , i = 1 + s p r_1 , 
, otherwise ; s = 0, ...,p — 2 . 



The R p r representation has a remarkable property 



R p CL 



"J^s p a — b [ I p b 



e.g. 



1Z \ I , 



(25) 



where TZ p a-b ( ) is formally coincided with R p a-b and has the entries of the identity matrix 



(26) 



I p b and zero matrix Q p b instead of 1 and 0. An action of R p r on the basic vectors read 

P-2 

R P - e p r-i( p _i) = -^e 1+spr -i , R pr ej = e j+1 , j = l,...,p r ~ 1 (p-l) - 1 , 

s=0 

The corresponding modular equations could be reduced 

rrij = q mod (n) , q = 0, n — 1 ; p m[ ? ' = p ^r-i^.^ = mod (n) 



2 



p = p m 



(y) 



i=i j=i 
j=p r - 1 {p-l)-l 



(<?) 

rrij = p 



- ... = p = mod (n) , 

i=i 

P r - x (p-i) f-^ip-t) 



E 



(?) 



i=p r - 1 (p-i)-2 



p m i = ^ m °d ( n ) 

i=P r " 1 (p-2)+2 



E(g) 
%- +p r-l (p _ 2 ) 



mod (n) . 



All equations are coincident if p = mod (n). Hence it follows N pT ^ p ~ 1 ^ = p. 

In Appendix |A] we give a straightforward solution for corresponding modular equations 
for the rotation axis Cg. 



D. Rotation symmetry k = Pi 1 p 2 2 , Pi ^ P2, r i, r 2 > !• 



This symmetry operation has an irreducible representation R p ri p r 2 through a unimodular 



matrix det R p T ^ p ^ = 1, dimi? p n p r 2 = p r ^ p T2 (p 1 — l)(p 2 — 1). Its p^ 1 p 7 2 2 - power is a 
diagonal block - matrix 



ri — 1 ro — 1 

Pi Pi 



ri — 1 ro— 1 



-R I) r l„ r 2 
Pi P2 



Pi P 2 



diag 



RpiP21 RpiP21 ■ ■ ■ 1 RpiP2 



which decomposes the action of lR p n p r 



. r-i — 1 ro — 1 
t p 2 2 



t = 1,2, ...,pip 2 , on the basic vectors 



ej, j = 1, 2, 1 ^2 2 ~ 1)G°2 — 1) into identical actions of ( R P1P2 ) on p^ 1 l p' 2 1 sets 



J"X — X m T2— 1 



of the basic vectors ej, j = 1, 2, (pi — l)(p 2 — !)• According to Section [Tl 1B| such action 



does not give rise to the nontrivial solutions of corresponding modular equations N, 



tfffc) 



where k = PiP 2 2 - 



E. Rotation symmetry k = Ylt^P^i r i — !• 



The generic case can be easily reduced to the previous Section [111D| and consequently 
does not give rise to the nontrivial solutions of corresponding modular equations. 



IV. CONCLUSION 

In the present paper we considered the colour lattice L d (n) with modular sublattices, 
when the only one crystallographic type of sublattices does exist and the only one of the 
colours occupies a sublattice, which is still invariant under /c-fold rotation C^. Such kind 
of colouring always preserves an equal fractions of the colours composed L d . The relation 
between the k - fold of the admitted rotation axis Ck and the number n of colours, which 
marked every modular sublattice, is given by the following formulae 



*(k) 



p , k = p r , r > 1 , 

1 > & = ]li=i 2 p? , n > l 



(27) 



which means that the n - colour lattice with modular sublattices allow to exist the crystal- 
lographic rotations for k = p r , r > 1 and n < p, where p is a prime number. A simple 
corollary of (p7|) gives an answer to the question: what is a minimal dimension d of the 
colour lattice L d (n) with n modular sublattices, which possess the rotation axis , k > n 



where p m j n is a minimal prime number which exceeds n. Both formulae ([27]), (p8|) extend the 
Hermann theorem of the crystallographic restrictions on the colour lattices with modular 
sublattices. 
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APPENDIX A: C 15 AND C 9 ROTATION AXES. 



First we illustrate the diagonalization (p0|) of irreducible representations As, Bi 5 for 
rotation axis Ci 5 by taking their successive powers. 



and 



.4 



15 



/ 4 — i? 5 

Ra —Rh 



A 



15 



3/ 



diag 



Rt 



Rt 



, t = 1,2,3,4,5 



B 



15 



V 



2 


2 


2 


-R 3 


R 3 


2 


2 


-R3 


2 


R 3 


2 


-R 3 


2 


02 


R 3 


-R 3 



B 



15 



5f 



diag 



R 



3 , 



R, 



R- 



R 



t = 1,2,3. 



Let us give a straightforward solution of modular equations for As irreducible representation 



.4 



15 



000 1 \ 
-1 1 

-1 1 
0-11 

-1 1 

1 -1 -1 1 
10-1 0-1 01 

\0 1 -1 -1 1 J 

Isometric transformation looks like 



As 


ei 


= e 6 , 


As 


e 5 = 


-e 2 


- e 6 


As 


e 2 


= e 7 , 


As 


e 6 = 


-e 3 


- e 7 


As 


e 3 


= e 8 , 


As 


e 7 = 


-e 4 


- e 8 


As 


e 4 


= ~ Ej=5 e i 


As 


e 8 = 




1 e i ' 



As-d 



q,n 



m 



(?) 



,(«) 



(?) 



e 2 + ™ ' - ml ' e 3 + m 



.(?) 



(?) 



(?) 



(?) 



M) 



(?) 



e 5 + ^8 + m l — m 



,(«) 



m 



(?) 



e 6 



(?) 1 (?) 
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Modular equations 
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are unsolvable for any modulus n. 

The corresponding representation Rg looks like 
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Rg ■ e 1 = e 2 , Rg ■ e 4 = e 5 , 
-R9 ' e 2 = ©3 , Rg ■ es = eg , 
i? 9 • e 3 = e 4 , i? 9 • e 6 = -e 1 - e 4 



It leads to the following modular equations 
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3m^ ; = 3 (m£" + <'j = 3 (m^ + m^) = mod (ra) , 
3 (mj° + mj + m^) = 3 (m^ + m? + m^) = mod (n) , 
3 (mf + m^ + mi 9) ) = 3 (m^ + rr$ + mjf^ = mod (n) , 



which have a solution iVf = 3. 
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FIG. 2. 



FIGURE CAPTIONS. 

FIG.l. Five - and two - colourings of the plane lattice with four - rotation axis C 4 . 

The first case n = 5 does not possess the permutation invariance of five sublattices, 
which have not had equal fractions in the unit cell. In opposite, the second case 
n = 2 does possess both these properties. 



FIG.2. 



Two - and three - colourings of the plane lattice with modular sublattices, which 
preserve the C3 rotation axis. Both lattices are built out of sublattices which 
preserved the full permutation invariance of the colours of equal fractions. 



